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ABSTRACT 



Through the use of asymptotic and perturbation methods, 
this thesis presents a theoretical study of the flow of a 
stratified fluid over variable topography as a model of the 
resonant interaction of a submarine's wake with a stratified 
fluid. Such resonant interactions may be able to produce 
significant upstream disturbances. The long time solution 
obtained in our model exhibits growth in time for the resonant 
case, indicating that perhaps nonlinear effects, balanced by 
dispersion, could cause significant upstream disturbances. 
These disturbances could allow the possibility of detecting a 
submarine by the internal waves its wake generates in the 
resonant case. 
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I. 



INTRODUCTION 



The nonacoustic detection of submarines is an area of 
intensive and well funded research. Although, it is not clear 
whether devices based on nonacoustic means of detection will 
become as useful as those based on acoustic means, the 
military importance of such devices demands pursuing such 
avenues of technological research. One possible nonacoustic 
means of detection is to attempt to detect the submarine 
optically through the use of a laser by examining the way the 
submarine reflects or absorbs blue-green light. Another more 
promising method is to try to detect the submarine by the heat 
it emits. However, an even more promising technique is the 
detection of submarines by the internal waves they generate. 

In the ocean a stable density stratification generally 
exists. As the submarine moves through the density- stratified 
ocean, it leaves behind a wake capable of generating internal 
waves. These internal waves result from the stabilizing 
effect of the stratification: as a water particle is 
displaced, it experiences a buoyant restoring force in the 
opposite direction from which it was displaced. 

Many theoretical and experimental studies have been 
conducted on the subject of the flow of a stratified fluid 
over a localized topography. For the two-dimensional flow of 
an inviscid, incompressible stratified fluid, most studies 
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have been confined to steady flow or linearized theory. 
Earlier studies concerned themselves only with a description 
of the downstream stationary lee -wave field. A question that 
has generated some controversy is the phenomenon of upstream 
influence. To what extent do disturbances generated by the 
topography through transient and nonlinear processes propagate 
upstream? 

McIntyre (1972) showed that the only significant upstream 
disturbances were weak second order long-wave motions 
generated by nonlinear interactions in the lee-wave. McIntyre 
(1972) extended the earlier theories of Benjamin (1970) and 
Ready (1971) by performing an expansion using as the small 
parameter the ratio of the height of the topography to the 
total fluid depth. Baines (1977,1979), on the other hand, in 
experimentally examining the flow of a continuously stratified 
fluid over topography, found upstream disturbances that were 
first order in the topography height. These upstream 
disturbances were not only generated by nonlinear processes 
over the topography but also travelled at the long- wave speed. 
Baines also observed that these upstream disturbances were 
very strong when the fluid approached resonance. Here, 
resonance is defined as the near coincidence of the basic flow 
speed and one of the free long-wave speeds. Grimshaw and 
Smyth (1986) presented a theoretical study of stratified fluid 
flow over localized topography for the two-dimensional case 
when the basic flow is near resonance. One interesting 
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feature of their results is the appearance of significant 
upstream disturbances as observed by Baines (1977,1979). 
Grimshaw and Smyth showed that solitary waves were formed when 
the localized forcing over the topography produced a positive 
disturbance upstream, and an oscillatory wavetrain was formed 
when the disturbance had negative polarity. 

Motivated by these results, this thesis presents a 
theoretical study of the flow of a stratified fluid over 
localized topography for the three dimensional case. 
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II. PROBLEM FORMULATION 



A. DIMENSIONAL PROBLEM 

We consider the three-dimensional flow of an inviscid, 
incompressible stratified fluid. Figure 1 (see page 5) 
provides a sketch of the coordinate system. We assume that 
the basic state has a constant horizontal velocity, V, from 
left to right with equilibrium pressure and density 
respectively as Po(z) and Pg{z) . Here and throughout this 
paper barred variables represent dimensional variables. We 
introduce f as the vertical particle displacement from some 
initial position at time t = 0 to the position at time t = T. 
Using the Lagrangian representation f(x,y,z,T) = Zf - the 
density, p(x,y,z,t), becomes ( z- f (x, y, z, t) ) . Therefore, the 
equations of motion are 



u- + V- + = 0 , 



( 2 . 1 ) 




( 2 . 2 ) 




(2.3) 



+ P? = 0 



(2.4) 
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Figure 1. The Coordinate System 
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w - 



El 

Dt' 



(2.5) 



where 



— = — ^ (V+u) -4 + v-4 + . (2.6) 

UE dT dx dy dz 

Here, u, v, and w are the velocity components relative to the 
basic state which has velocity (V,0,0) with respect to a 
right-handed Cartesian coordinate system. 

The bottom topography is given by z = f(x,y,t), where 
initially f(x,y,0) = 0. For large times the topography 

attains a steady value g(x,y) and hence we assume that 

f{x,y,t) - g(x,y) as t -» <». (2.7) 

We also assume that at all times the topography is localized 
where f (x,y,t) 0 as x^+y^ <». The Icinematic bottom boundary 

condition is then 

7 = f at z = f . (2.8) 
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We are assuming that this "topography" is created by the 
submarine's wake, which when viewed from the moving submarine, 
actually has a steady shape. In this thesis, we investigate 
the interaction of the wake with internal waves. Kallen 
(1987) has investigated this type of interaction in a linear, 
steady state model and comments in detail on the validity of 
assuming that the wake acts as rigid topography. Furthermore, 
from experimental results, Lin and Pao (1979) show the mean 
flow occurs at the wake boundaries implying that part of the 
flow travels around rather than through the wake. It is this 
flow around the wake that generates the internal waves. Our 
model, in fact, does not depend critically on the assumed form 
of "topography" the submarine creates and our purpose here is 
to investigate a time dependent model. 

The upper boundary, assumed to be a free surface, is 
located at z = H + i;(x,y,t), where H is the vertical distance 
from the surface to the submarine and r/(x,y,t) is the free 
surface displacement. One must remember to include not only 
the kinematic boundary condition at the upper boundary but the 
dynamic boundary condition as well. The upper boundary 
conditions take the form: 



f = at z = H + 7} , 



(2.9a) 



P^{H+7}) + Pp{H+rj) = p„{H) at z = H + r} . 



(2.9b) 



7 



Here, Pj, is the perturbation to the pressure p,, of the basic 
state . 

B. DIMENSIONAL SCALING 

The solution of practical problems in fluid mechanics 
requires both theoretical developments and experimental 
results. By grouping significant dimensional quantities into 
dimensionless parameters it is possible to present results in 
a compact form which is applicable to all similar situations. 

In the problem, the characteristic horizontal length scale 
L is determined by the ratio U/N, where U is a characteristic 
horizontal velocity and N is the Brunt -Vaisala frequency at 
the location of the submarine. We shall take the 
characteristic horizontal velocity U to be given by N,H, (U = 
N,H) , where N, is the maximum Brunt -Vaisala frequency and H, 
the characteristic vertical length, is the depth of the 
submarine. Typically Nj is 10'^ sec'' and H is 200 meters (m) , 
so U is 2 m/sec. This is a typical order of magnitude 
patrolling speed for a submarine at that depth. Of course, 
0(1) variations in N, and H can be made without changing our 
analysis . 

From the mass conservation equation and with the imposed 
horizontal length and velocity scales, we see that the 
vertical velocity scale, W, is UH/L. The horizontal length 
and velocity scales are also used in determining the 
characteristic time scale, T, where T is defined by T = L/U = 



8 



1/N. Lastly, the characteristic pressure scale is defined as 
p = PigH where p, is a characteristic density. These scales 
define the parameter = N,^H/g, which is a dimensionless 
measure of the ratio of the buoyancy force to the 
gravitational force. 

C . NOND IMENS I ONAL I Z AT ION 

By using the various dimensional scales introduced in 
section B, we shall nondimensionalize the equations of motion 
(2. 1-2. 6) along with the boundary conditions (2. 8-2. 9b). 
However, before deriving the nondimens ional equations of 
motion and boundary conditions, we will consider several 
related dimensional equations that will assist in simplifying 
the nondimens ional derivation of equations (2. 1-2. 6) and 
boundary conditions (2. 8-2. 9b). We first introduce the 
following nondimens ional variables: 



_ Po 


z = ^ 




H 


N 


r - i . 

^ H 


■ IT, ' 




( 2 . 10 ) 



To nondimensionalize the hydrostatic equation, we begin with 
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( 2 . 11 ) 



dPo 

dz 



= ~9Po- 



Using (2.10) yields 



P'o (z) = -Po- 



( 2 . 12 ) 



The dimensional Brunt -Vaisala frequency is defined as 



N'^ 



g <^Po 

To ^ 



(2.13) 



where dp^/dz < 0. Using (2.10), we obtain 



nIh 

~9~ 






1 dPo 
Po dz 



(2.14) 



but /3 = N,^H/g, therefore (2.14) simplifies to 



= - 



dPo 
dz ' 



(2.15) 
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In nondimensionalizing the Lagrangian representation of the 
particle's vertical displacement, we start with 



z 



o 



z - X ix,y, z, t) . 



(2.16) 



Using (2.10) yields 



z 



O 



z - f {x,y, z, t) . 



(2.17) 



The last equation needed before considering the equations of 
motion and boundary conditions is the dimensional pressure 
equation : 



where Pp = p,gHi8q (x, y , z , t ) and /3 is defined after (2.14). 
Thus, q(x,y,z,t) is the dimensionless perturbation to the 
hydrostatic pressure p„. 

Before deriving the nondimens ional equations of motion and 
their associated boundary conditions, we first introduce the 
following nondimens ional variables: 



p(x,y,z,t) = p^gHp^iz) + p^gH^q(x, y, z , t) 



(2.18) 
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X = 


X 

L ' 


y = 


y 

X ' 


z = -5 , 

H 


u = 


u 

U ' 


V = 


V 

U ' 


V = 


V 

u ' 


w 

W 


Po = 


Po 

Pi 


t = 


t 

~f ' 


f = 


7 

H ' 


II 

1 


f = 


1 

h ' 



where L, H, U, W, p,, and T are the characteristic dimensional 
scales defined in section B and h is a characteristic 
amplitude of the topography. 

We begin with the mass conservation equation (2.1) . Using 
(2.19) and remembering that W = UH/L, (2.1) becomes 

+ w, = 0 . (2.20) 

To assist in deriving the nondimens ional form of the x- 
component of the momentum equations, given by (2.2), we shall 
now derive the nondimensional form for the convective 
derivative D/Dt: 



— = — ^ (V+u) + v_^ + w_^. (2.21) 

DT dT dx dy dz 



Using (2.19), (2.21) becomes 
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(2.22) 



D 

'ut 



1 a 

T TE 



Li 







d 

■?i 



but T = L/U, therefore (2.22) simplifies to 



_D ^ 1 D 
Dt T Dt' 



where 



D 

~Dt 




(V+u) 



d 



+ V 






+ 




Substituting (2.17), (2.18), and (2.23) into (2.2) 

(2.19), we obtain 



PxPo^^ 



n RPH 

^ ~T Dt 



Px9W 



= 0 or. 



Po 



(z- 



■ f)£H 

^ Dt 



+ g, = 0. 



A similar analysis for the y- component of the 
equations (2.3) yields 



(2.23) 



(2.24) 



and using 



(2.25) 



momentum 
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(2.26) 



( „ Dv 






To derive the nondimens ional equation for the z- component of 
the momentum equations (2.4), we begin with 



+ pJz-09^ Pe 



= 0 



(2.27) 



Substituting (2.17), (2.18), and (2.23) into (2.27) and using 

(2.19), we obtain 



Po(z-f) + Po^^-0 9 + 9Poz + 9^9, = 0 . 



(2.28) 



Using T = L/N,H and W = N,hVl, (2.28) simplifies to 




iz-n — 

^ Dt 



+ Po(z-09 + srp,^ 



+ 90q, 



0 . 



(2.29) 



Using (2.12) and the definition of /3, )S = N,^H/g, (2.29) 

simplifies to 
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(2.30) 



e^pjz-^)^ + pjz)) + g, 



= 0 , 



where 



e 



H 

L ■ 



(2.31) 



Lastly, we derive the nondimens ional vertical velocity 
from its definition (2.5). Using (2.19) yields 



w = 



Rl 

Dt ■ 



(2.32) 



We now consider the kinematic bottom boundary condition 



f = f at z = f. (2.33) 



Using (2.19), we obtain 



f = otf at z = Off, 



(2.34) 



where 
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a 



h 

H 



(2.35) 



For the free surface kinematic boundary condition, 



f = T) at z = H + 7} . (2.36) 

We use (2.19) to obtain 

f = Tj at z = 1 + t}. (2.37) 

Next, we nondimensionalize the free surface dynamic boundary 
condition in (2.9b): 

p^{H+^) +p^(H+tj) =p^{H) at 'z = H + ^. (2.38) 

Using (2.19) and (2.18), (2.38) becomes 

p^(l + ij)+j8g = p„(l) at z = l + tj. (2.39) 

To recapitulate our nondimensional versions of equations 
(2. 1-2. 6) and the surface boundary conditions (2. 8-2. 9b), the 
nondimensional equations of motion are 
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Du 

Dt 



Dv 



+ <Iz = 0, 



where e = H/L and 



D _ d 
Dt Jt 



with boundary conditions 



p„(l+T/) + jSg = p^(l) at z = 1 + 7), 



+ = 0, 


(2 . 40) 


+ 

iP 

II 

O 


(2.41) 


o 

II 

-h 


(2.42) 


O 

II 

-1- 


(2.43) 


Dt 


(2.44) 


+ vj- + wj-. 


(2.45: 


z = af, 


(2.46) 


z = 1 + TJ, 


(2.47a) 


Z = 1 + TJ, 


(2.47b) 



where a = h/H. 

An exact analytical solution of these nonlinear partial 
differential equations is not available. However, we shall 
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obtain approximate solutions for the case when e defined by 
(2.31) and a defined by (2.35) are both small positive 
constants. Since s = H/L = N/N, , the condition e << 1 is 



satisfied when 


the density changes 


rapidly 


across a 


thin 


thermocline 


• 


However, there is 


some ambiguity in 


the 


definition 


of 


the thermocline. 


We will 


consider 


the 


thermocline 


to 


be the region where 


density 


variations 


are 



determined chiefly by temperature variations. This means that 
the submarine is usually submerged to a depth below the 
thermocline. Is this a valid assumption? To avoid acoustic 
detection the submarine submerges to depths below the mixed 
layer or layer depth. Depending on the thickness of the 
thermocline the submarine can be either in the thermocline or 
below it. We shall consider only the case for which the 
submarine is below the thermocline in which the parameter a, 
the ratio of the characteristic amplitude of the topography to 
the submarine's depth defined by (2.35) , is also small in this 
situation. 
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III. THE ASYMPTOTIC EXPANSION OF THE DEPENDENT VARIABLES 



To study the behavior of the internal waves excited by a 
small amplitude disturbance in the presence of a uniform 
oncoming flow, as when a submarine moves through the water at 
constant velocity V, we perform an asymptotic expansion on the 
variables q, u, v, and tj in terms of the small parameter a. 
This is reasonable since the small parameter a enters the 
bottom boundary condition in an analytic way and since the 
solutions arise by analytic processes, it is natural to expect 
the solution to exhibit an analytic dependence on the 
parameter a. Since the disturbance is 0 {q!) , we shall 
initially suppose that the response is also 0(o;) . Therefore, 
we introduce the following asymptotic expansions 

f = oi^Jx.y, z, t) + + . . . 

g = oiq„{x,y,z, t) + a^g, + . . . 

u - aUg{x,y, z, t) + oi^u^ + . . . (3.1) 

V = av^ (x, y, z, t) + a^v^ + . . . 

rj = oit)^{x,y, z, t) + + . . . 

in which we consider the limit process as a ^ 0. 

Substituting (3.1) into the continuity equation (2.40) , we 
obtain 
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= 0 . 



(3.2) 



a 

~53t 



[auj 



-I- 



a 

■5y 



[avj 




D 

Dt 



(ar.)] 



+ O(cx^) + . . . 



Dividing by o' and taking the limit as of 0 , (3.2) 

to 






V + 
^ay 



= 0/ 



where 



r 

Dot 



and 



Do 

Dot 



a 

TE 







To show that Do/D^t is the 0(1) part of the 
derivative given in (2.45), we collect only the 0( 
the expansion for D/Dt given in (3.6) below: 



D 

Dt 



a 

TE 






+ oiv„ 



■3y 



0{oi^) 



simplifies 

(3.3) 



(3.4) 



(3.5) 

convective 
1) terms in 



(3.6) 
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Substituting (3.1) into the x- component of the momentum 
equations (2.41), we obtain 



p„{z-^) (au„) + aq^ + 0{a^) + . . . = 0 . (3.7) 



Expanding p„(z-^) in a Taylor Series and using (3.1) yields 



Pjz-r) = Poiz) - p'jz) [(xU + O(a^) + (3.8) 



Substituting (3.8) into (3.7), we obtain 



<xp„(z) + aq^ + O(oi^) + . . . = 0. (3.9) 

^ rk ^ 



Dividing by a and taking the limit as a ^ 0, (3.9) becomes 



Po(2) 



DX 



<?<„ = 0- 



(3.10) 



A similar analysis for the y- component of the momentum 
equations (2.42) yields 
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(3.11) 



D V 

p (z) ” ^ + a =0. 

' n <- ^oy 

O ^ 

For the z-component of the momentum equations (2.43), we 
substitute (3.1) into (2.43) to obtain 

0 LeP--l-l + .i [p^(z-f) - p„{z) ] + org^ + O(a^) + . . . = 0. (3 .12) 

D^t P 

Expanding p^(z-^) in a Taylor Series as in (3.8) and 
substituting into (3.12), we obtain 

- -^P^(Z) + Oiqoz + + . . . = 0. (3.13) 

D^t P 

Dividing by of and taking the limit as of ^ 0, (3.13) simplifies 

to 



-^p^(z)r. + g,. = 0. (3.14) 

Notice that the term has been neglected. As previously 
discussed, in order for our theory to work we require e << 1. 
Therefore, we see the term is small. The fact that the 
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term is small by itself is not justification in dropping this 



term from (3.13) . In order to justify the dropping of the 
term, we shall assume that or£^ = o(a) as a ^ 0, or that = 
o(l) as O' -* 0 . Now as one ta)ces the limit or -* 0 in (3.13) , we 
clearly see that -* 0. In continuing our derivation of the 
z- component of the momentum equations, we substitute (2.15) in 
for dpo/dz, so that (3.14) becomes 



Next, we consider the surface boundary conditions. We 
begin by substituting (3.1) into the Icinematic bottom boundary 
condition (2.46) to obtain 



Note that the 0(o;) bottom boundary condition is applied at z 
= 0. Substituting (3.1) into the free surface kinematic 
boundary condition (2.47a) results in 




(3.15) 



Q!^ = cuf + O(oi^) + . . . at z = af. 

o 



(3.16) 



Taking the limit as or ^ 0, we find 



= f at z = 0. 



(3.17) 
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= 017) „ + O(oi^) + , , , at z = 1 + 017)^. (3.18) 

Taking the limit as a ^ 0, (3.18) becomes 

= rj, at z = 1. (3.19) 

Again the 0(o!) free surface condition is applied at z = 1. 

The last boundary condition to consider is the free 
surface dynamic boundary condition (2.47b) . Substituting 
(3.1) into (2.47b), we obtain 

p„(l + rj) + Q!/8g„ + O(a^) + . . . = p^(l) . (3.20) 

Expanding p„(1 + t;) in a Taylor Series and using (3.1) gives 

P„(l + T/) = p^(l) + Otp'„{z) 7)^ -<r 0{(x^) + — (3.21) 

Substituting (3.21) into (3.20), we find 

p„{l) + Oipl iDVo + + O (q;^) + • . • = Po ( 1 ) . (3.22) 
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Canceling the Po(l) terms, dividing by a, and taking the limit 
as q; -» 0, (3.22) becomes 



Po (z) »?„ + = 0 at z = 1 



(3.23) 



but using (2.12) simplifies the free surface dynamic boundary 
condition to 



= P^Vo at z = 1. (3.24) 

In summary, we see from substituting the asymptotic 
expansion (3.1) into (2.40-2.45), (2.46), and (2.47a,b), that 

the 0(q;) problem is: 



= 0^ (3.25) 

+ g«. = 0, (3.26) 

Po^^)^ ^ <Joy = 0, (3.27) 

g« + = 0, (3.28) 
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(3.29) 






where 



Do 

Dot 



d 




(3.30) 



with boundary conditions 



II 


at 


z = 0, 


(3.31) 


L = -no 


at 


z = 1, 


(3.32a) 


= PoVo 


at 


z = 1. 


(3.32b) 



Combining (3.32a) and (3.32b) gives the single boundary 
condition at the free surface: 



&<io = PcL -at z - 1. 



(3.33) 
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IV. DEVELOPMENT OF THE TWO-DIMENSIONAL WAVE EQUATION 



A. EIGENVALUE PROBLEM 

In order to solve (3.25-28), (3.31), and (3.33), we look 

for solutions of the form 



u, = 4>' (z) 

0' (z) 

= DeHKX*LY-OT) 

ct>' {z) 



(4.1) 



where k = (K,L) is the horizontal wave vector, Q is the 
frequency, and c is the long-wave phase speed relative to the 
basic state at rest. Substituting (4.1) into (3.25), (3.26), 
and (3.27) and solving the system of linear equations, we find 
the dispersion relationship to be 



Q = c I K I where | k | ?^ 0 . (4.2) 

Furthermore, substituting (4.1) into (3.28), (3.31), and 

(3.33) results in 
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c 



(4 . 3a) 



with boundary conditions 



(f>(z) = 0 at z = 0 , 



(4.3b) 



0(z) = {z) at z = 1. 



(4,3c) 



The equations (4.3a-c) formally pose an eigenvalue problem 
where c'^ is the eigenvalue and 4> the eigenfunction. Note that 
the dispersion relationship in (4.2) was obtained for wave 
motion in a stratified fluid at rest where V = 0 and D^/D^t = 
d/dt. In the problem we are considering, the dispersion 
relationship would differ from that found in (4.2) but the 
eigenvalue problem would remain the same. Hence, our purpose 
here is to show the existence and character of the solutions 
that define the eigenvalue problem. The solutions have the 
horizontal character of plane waves travelling in the (x,y) 
plane with a vertical structure defined by the eigenfunctions, 
0(z). We shall see that the general solution to (3.25-3.28), 
(3.31), (3.33) can be obtained in terms of a superposition of 
normal modes whose depth dependence is determined by the 
</>s(z), the eigenfunctions of the following eigenvalue problem: 
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(4.4a) 



(f)^{z) = 0 at z = 0, 



(4.4b) 



4>,{z) = I3c^(f)^^(z) at z = 1. 



(4.4c) 



Here is the eigenvalue, where is the long- wave phase 
speed relative to the basic state at rest. 

B . ORTHOGONALITY 

The eigenfunctions <t>s(z) satisfy the orthogonality 
condition 



where is defined by the left side of (4.5) for r = s. To 
prove (4.5) , we let and X, be eigenvalues with corresponding 
eigenfunctions and 0, where X^ = c/^ and X, = c,'^. The 




(4.5) 



differential equations satisfied by the eigenfunctions are 






(4 . 6a) 



^PoK'>' + KPo^^4>r = 0 - 



(4.6b) 
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In addition, both </>s and 0, satisfy the same set of homogeneous 
boundary conditions defined in (4.4b,c). Multiplying (4.6a) 
by and (4.6b) by yields 



K^Po<i>sz)'4>r + KKp„N^4>s4>r = 0, (4.7a) 

K (PoK )'<t>s ^ KK PoN^4>s <t>r= 0. (4.7b) 

Subtracting (4.7b) from (4.7a) results in 

KiPoK)'^r - KiPo<Pn)'^s - 0. (4.8) 

Through application of the product rule 

K(Po<t>s^)<t>r = KUPcK^r)^ - PoKK^> ( 4 . 9 ) 

(4.8) becomes 

K^Po<i>sz<i>r)' - KPo<t>sz<t>n ~ K ^ Po4>n <t>s ^ ' + \PoK<t>sz = 0 - (4.10) 

Integrating (4.10) yields 
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rl ^ ^ 

- V) Po4>szK^^ = I “ KPo^sz I • 

JO 0 0 



(4.11) 



Evaluating the right hand side of (4.11) and remembering 0^(1) 
= |SX;' 03 ,( 1 ), </>,(!) = ^X/‘03,(l), and <t>,{0) = 0,(0) = 0, (4.11) 

becomes 



(X, - X,) 




dz = 0 . 



(4.12) 



If X, it follows that 

Jq pAj 0;, dz = 0 . (4.13) 

In other words, the z derivatives of the eigenfunctions (0, 
and 0,) corresponding to different eigenvalues (X, and X,) are 
orthogonal with weight Therefore, the normal modes do 

satisfy the orthogonality condition in (4.5). 

C. TWO-DIMENSIONAL WAVE EQUATION 

Since the eigenfunctions, 0,(z), form a complete set, we 
can represent u^, v„, i*o, and qo by a generalized Fourier series 
expansion. From (4.1) we know the general form of the 
solutions and hence we put 
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(4 . 14a) 



To derive a 



We take the 
by PoCs^0rz(z) 
obtain 



Integrating 



Co 






= 



Po ( ^SZ < ' 



0 



Uo = -Y^D^{x,y, t)^^^(z) , 
0 

oo 



relationship between and B^, we begin 



Co = t) 4 >s( 2 ) . 

0 



derivative of (4.15) with respect to z, 
and use the orthogonality condition in 



^sAs 



= c. 



! f P oC OZ 4* SZ ■ 

J 0 



(4.16) by parts then gives 



(4 . 14b) 
(4 . 14c) 
(4 . 14d) 

with 

(4.15) 

multiply 
(4.5) to 



(4.16) 
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(4.17) 






cl?o^o^sz\ - 

0 *'0 



Using (4.4a), (4.17) becomes 



= 



-sPo^szl * 

O Jo 



(4 . 18) 



Now, we multiply (4.14b) by (/)„ and from the orthogonality 
condition (4.5), we obtain 









,dz. 



Integrating (4.19) by parts, we obtain 



0 Jo 



Using (3.28), (4.20) becomes 



= ^o 4 »si + pPo^o^^<|)sC(z. 
0 Jo 



Combining (4.18) and (4.21) yields 



(4.19) 



(4.20) 



(4.21) 
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(4.22) 



cl?o^o^sz\ + ^sBs - <7o<l>J = 



Substituting (3.33) and (4.4c) forq^d) and 0^(1) respectively 
and remembering that f„(0) = f(x,y,t), (4.22) simplifies to 



Bs = F„ 



(4.23) 



where 



= c|po(0)<t>5^(0) f (x,y, t) . (4.24) 

Note that Fs(x,y,t) is proportional to the bottom topography 
f (x,y, t) . 

Next, we substitute (4.l4a,c,d) into the continuity 
equation (3.25) to obtain 




dx 







D t 

0 ■^o‘' 



= 0 . 



(4.25) 



Multiplying (4.25) by and applying the orthogonality 

condition (4.5) then gives 
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(4.26) 






D^t dx dy 



where Dp/D^,t is given in (3.30) . 

To derive the two-dimensional wave equation, we use the x 
and y components of the momentum equations (3.26) and (3.27) . 
Substituting (4.14b,c) and (4.14b,d) into (3.26) and (3.27) 
respectively, we obtain 



Multiplying (4.27a) and (4.27b) by and applying the 

orthogonality condition (4.5) then yields 




(4.27a) 




{A. 27b) 




(4.28a) 




(4.28b) 
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We next take the x derivative of (4.28a) and the y derivative 
of (4.28b) and add the two results to obtain 



D„t 



dD, dE^ 
T3E 






+ = 0 



(4.29) 



Substituting (4.23) and (4.26) into (4.29) finally yields 











(4.30) 



where 



and 





(4.31) 




a 

■?t 




(4.32) 



Equation (4.30) is the inhomogeneous two-dimensional wave 
equation whose solution will be obtained in Chapter V. 

To summarize, we see that substituting (4.14a-d) into 
(3.25-3.27) implies 



36 



B, 



(4.33a) 



= A. 






D^t bx dy ' 



(4.33b) 






= CjVM A, + F ] , 



(4.33c) 



where 



I,F, = c^pJ0)4>,A0) fix.y, t) , (4.33d) 



and 



Do 

Dot 



d 

TE 



+ V 



d 

■53E' 



. 



(4.33e) 



(4.33F) 



Here, since f vanishes at t = 0, Fj(x,y,0) = 0 and from (2.7) 
we have 



F^(x,y, t) ~ G^(x,y) 



as t -* 00 , 



(4 . 34a) 



where 
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(4.34jb) 



I,G, = cjpJ0)4>s,{0)g{x.y) . 



The initial conditions for (4.14a-d) are that A^, B^, D^, and 
Ej all vanish at t = 0 . 
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V. SOLUTION TO THE TWO-DIMENSIONAL WAVE EQUATION 



A. METHOD OF DESCENT 

To solve (4.33c) we first introduce the following change 
of variables: 



^ = X - Vt, 

y = y, (5 . la) 

a = t. 



Derivatives with respect to x, y, t can be expressed as 
derivatives in the new variables y, a as follows: 



d 

~53c 

d 

■5y 

d 

Tt 



d 

d 

■5y' 

d 




(5.1b) 



Substituting (5.1a,b) 



into (4.33c) yields 



d% 2 r 

— ^ - c; [ f + — i 

da^ dy^ 






(5.2) 
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where 



.2 r , a^ 



o) = c, [ -_ + __ ] F,(^+Vo,y, a) . 

a r ay^ 



(5.3: 



We shall solve (5.2) by beginning with the solution to the 
following three-dimensional wave equation: 



a^u ^ ^2 [• a^u , a^u . a^u 



at^ 



ax^ ay^ az^ 



] , 



u{x,y,z,0) =^{x,y,z), 



u,(x,y,z,0) ='^{x,y,z). 



(5.4a) 

(5.4b) 

(5.4c) 



Using Kirchoff's formula (see Strauss 1992), the solution to 
(5 . 4a- c) is 



u{x^,y„, z^, 



47rc^t, 



- T^(x,y, z)dS 

o 



a 



47TC 



i_ r r 4> (x,y, z) dS 
c't, ;j 



(5.5) 
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where 



D = [ {x,y,z) : (x-x^) ^ + (y-y„) ^ ] . (5.6) 



Note that the value of u (x„, y„, z„, t„) depends only on the values 
of '^(x,y,z) and <f>(x,y,z) on the spherical surface D. 

By the method of descent, we now obtain the solution to 
the following two-dimensional wave equation: 



^ ^ a^u ^ 


(5.7a) 


de dx^ dy^ 




u (x,y, 0) = <l> (x,y) , 


(5.7b) 


u,(x,y, 0) = '^^(x,y) . 


(5.7c) 



We regard the solution u(x,y,t) of (5.7a-c) as a solution to 
the three-dimensional problem (5.4a-c) with initial data 
independent of z. By Kirchoff's formula (5.7), we have 



u(0,0,t) = — ^1— r \'ir{x,y)dS 

47TC^t IJ 



(5.8) 



TE 



47TC^t . 



4>(x,y) dS 



D. 
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where 



D, = [ {x,y,z) : ] . (5.9) 



This is the solution to (5.7a-c) but it can be simplified. To 
simplify (5.8), we must express the surface integral as an 
ordinary double integral. It is twice the integral over the 
top hemisphere z = (c^t^-x^-y^) where (x,y) ranges over the 
region D, = { (x,y) : x^ + y^ s c^t^ }. By definition 



dS = I X I dxdy, 



(5.10a) 



where 




(5.10b) 



and 



\R.^Ry\ = , .. i: 

/(c^t^-x^-y^) 



(5 . 10c) 



Substituting (5.l0a-c) into (5.8) yields 
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u(0, 0, t) 



'^(x,y) 



+ 







dxdv 

\l c^t^-x^-y^ 



dxdy 



(5.11) 



This is the solution at the point (0,0,t) . Therefore, the 
solution at a general point is 



Note that (5.12) shows that the value u(x^,yo,to) depends on 
the values of <i>(x,y) and '^^(x,y) inside the disc defined in 
(5.13). This illustrates the well known fact that Huygens' 
principle fails in two dimensions. We shall use (5.12) in the 
solution to the nonhomogeneous two-dimensional wave equation 
discussed in the next section. 




(5.12) 




where 



D, = [ (x,y) : (x-x^)^ + (y-y„) ^ s 



(5.13) 
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B . DUHAMEL ' S PRINCIPLE 



We now consider the following nonhomogeneous two- 
dimensional wave equation 



d^u _ ^2 [ ^ 

dt^ dx^ dy^ 



g(x,y, t) , 



(5.14a) 



with homogeneous initial conditions 



u (x,y, 0) =0, 


(5.14b) 


u,(x,y, 0) ^0. 


(5.14c) 



By Duhamel ' s principle, the solution to (5.14a-c) is given by 



u(x,y,t) = C7(x,y, t, s) ds 



(5.15) 



where for each s a 0, U(x,y,t,s) is the solution of 



^ - cM ] 

b bx^ by^ 



0 for tz s, 



(5.16a) 



with initial data prescribed on the plane t = s such that 
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U(x,y, s, s) = 0, 



(5 . 16b) 



U,(x,y,s,s) =g{x,y,s). (5.16c) 

Since the wave operator is invariant under translations in 
time t (see Strauss 1992) , we can shift the time t by a 
constant and still have a solution. We let V(x,y,t,s) 
U(x, y, t+s, s) , so that V(x,y,t,s) is a solution to 



K + V^) = 0 for t a 0, (5.17a) 

with initial data prescribed on the plane t = 0 such that 

V(x,y,0,s) ^ U(x,y,s,s) =0, (5.17b) 

V,(x,y,0,s) =U,{x,y,s,s) =gr(x,y,s). (5.17c) 

Using Kirchoff s formula (5.12) to obtain a solution to 
(5.17a-c) then yields 



V(x,y, t,s) = 



27TC 



J 



g(x,y, s) 



-Id, y/ c^t^- (x-x) (y-y) 



dx dy 



(5.18) 
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where D, = { (x,y) : (x-x)^ + (y-y)^ s }. Since our purpose 

here is Co solve the nonhomogeneous two-dimensional wave 
equation in (5.14a-c) by Duhamel's principle (5.15), we can 
again shift the time t by a constant so that U(x,y,t,s) = 
V (x, y , t - s , s ) . This last expression for U(x,y,t,s) together 
with (5.18) then implies 

U(x,y,t,s) = r r g(x^y, s) ^ ( 5 . 19 ) 

"Id, s/ c^ ( t-s) (x-x) (y-y) ^ 



where 



D, = [ (x,y) : (x-x)^ + iy-y)^ s cMt-s)^ ] . (5.20) 

Now using Duhamel's principle (5.15), we obtain 



u(x,y,t) = 



g(x,y, s) 



= = 

•7)/ V ( t-s) (x-x) {y-y) ^ 



cZx dy 



ds (5.21) 



the solution to (5.14a-c) where D, is defined in (5.20). 

In deriving the solution to the nonhomogeneous two- 
dimensional wave equation (5 . 14a- c) , we have solved (5.2). We 
change the names of the dependent and independent variables in 
(5.21) to obtain the solution of (5.2): 
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fJx+Vs,y, s) 



^ 

D. \J c^ (a-s)^- ix-^)^- iy-y)'^ 



dx dy 



ds 



(5.22: 



where D„ = { (x,y) : (x-f)^ + (y-y)^ s c^(a-s)^ }. Letting f = 

X - Vt where a = t, (5.22) becomes 



Ajx,y,t) = 






f^{x+Vs,y, s] 



(5.23) 



( t-s) (x- (x-vt) ) (y-y) ^ 



dxdy 



ds 



where Dj = { (x,y) : (x-(x-vt))^ + (y-y)^ s Cs^(t-s)^ }. 

To derive the solution for the nonhomogene ous two- 
dimensional wave equation, we began by solving the homogeneous 
three-dimensional wave equation using Kirchoff's formula 
(5.5). Next, we used the method of descent to solve the 
homogeneous two-dimensional wave equation by reducing 
Kirchoff's surface integral formula to an ordinary double 
integral (5.12) . We developed the solution to the homogeneous 
two-dimensional wave equation in order to solve the 
nonhomogeneous two-dimensional wave equation by applying 
Duhamel ' s principle (5.21). This gave the solution of the 
wave equation in (5.2) in the form given above in (5.23) . In 
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Chapter VI we will examine the solution to (5.23) for long 
time . 
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VI. LONG TIME SOLUTION 



A. C/ 

Before looking at the behavior of (5.23) as t -> c» (where 
the topography f 3 (x+Vs,y,s) has a general form), we shall 
first introduce a specific moving point source for the 
topography, so that 

f^{x,y,s) = 6 {x-x^) 6 {y-yj g(s) . 

Consequently 

f^{x+Vs,y, s) = 6 ix+Vs-x^) 6 {y-yj g(s) , 

where g is a non-negative function satisfying 

g(0) =0 and g(s) = l for s ^ s^. (6.1c) 

If we can solve (5.23) by substituting (6.1b) in for 
fj (x+Vs, y , s) , then we can solve (5.23) for any form of 
topography by applying the principle of superposition. 

To solve (5.23), we begin by substituting (6.1b) into 
(5.23) to obtain 



(6.1a) 



(6.1b) 
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A^{x,y, t) 






5 {x+Vs-x^) 6 {y-yj g{s) 

W [x-{x-vt) ]^-(y-y)^ 



( 6 . 2 ) 



dx dy 



ds. 



We let X = (x+Vs-x^,) and Y = iy-Yo) and substitute into (6.2) 
to obtain 



AJx,y,t) = 



1 

^ r 


r 


6 {X) 6 {Y) g{s) dXdY 


27Tc4 


o; 


\jcl ( t-s) [x+v( t-s) - {x-xj ] ^- [ Y- (y-y„) ] ^ 



(6.3) 



where 



D, = 



{X, Y) 



[X’+V( t-s) - (x-x„) ] ^ + [Y+y^-y]^ s c^^(t-s)^ . 



(6.4) 



In order to solve (6.3) we must determine those s for which 
(6.4) is satisfied. Letting t = (t-s) for (t-s) a 0, (6.4) 
simplifies to 



(c,^-V^)t^ + 2V{x-x„)t - [(x-xj^ + iy-yj^] a 0. (6.5) 
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Note that X and Y are not included in (6.5) since the 
integration in (6.3) with respect to X and Y makes both X and 
Y become zero. Solving (6.5) for t yields 



= 



-V{x-x) 






C. (X-X„ 



(c,-u^) (y-yj 






:6.6) 



The values of and t. will depend on whether c^^ > or c^^ < 
V^. We shall examine both cases in the following sections. 

1. C,^ > 

When the long -wave phase speed is greater than the 
velocity of the submarine, we see from (6.6) that t+ > 0 and 
T. < 0 . Factoring (6.5) yields 



(c^-V^) {t-tJ (t-tJ aO (6.7) 

where and t. are defined in (6.6) . Since the left hand side 
of (6.7) must be greater than or equal to zero, r must be 
greater than or equal to (t a t+) . Hence, s must be less 
than or equal to (t-r^.) . Therefore, the upper limit of s in 
(6.3) becomes (t-r+). Now, integrating (6.3) with respect to 
X and Y yields 
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A^{x,y,t) = 



2ttcJc;-V^} 



1/2 



i:- 



g{s) ds 



^(( t-Tj -s)((t-Tj -s) 



( 6 . 8 ) 



where and t. are defined in (6.6). Using (6.1c), we break 
(6.8) into the following parts 



A^(x,y, t) = J, + Jj, 



(6 . 9a) 



where 



and 



1 P. gis) ds 

1 r' ds 

2TTCjC^-V^)''^k -s)((t-T.) -s) 



(6.9b) 



(6.9c) 



We first look at the behavior of I, as t -* oo. We see 
from (6.9b) that for long time, I, is 0(l/t). Therefore, we 
can neglect this portion of Aj.(x,y,t) . Now in evaluating 
for long time, we let t = (t-s) so that Ij simplifies to 
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^2 = 



27TC, ( C, -V^) 



2 \ 1/2 



i: 



dr 






(6.10: 



We complete the square of the denominator of (6.10) to obtain 



^2 = 



2itcAc:-v 



2 \ 1/2 



r- 



dr 



\/ ( r-a) 



(6.11: 



Here a = (t + + t.)/2, and = (t + -t.)^/4. Letting bu = (r-a) and 
substituting into (6.11) then yields 






27TC (c,-V^) 



1/2 



I 



du 



v/u'-l 



(6.12: 



Integrating Ij, we find 
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1 


,n| . 


y / 2 

(t-s„-a)2 , I 




’ Jb 


J J 



2Trc,(c^-V^) 





^ >1 


1/2 


ml * 


(T,-a)^ _ 1 


1 


' b 


b^ 


1 



(6.13) 



For long time t, Ij is 0(ln t) . 

In the case where > V^, we see that the long time 
solution for A,(x,y,t) is 0(ln t) . Next we examine the case 
for which < V^. 

2 . 

When the long -wave phase speed is less than the 
velocity of the submarine, we see from (6.5) that x must be 
greater that x„ since t is greater than or equal to zero. We 
must reexamine (6.6) to see if a positive real value for t is 
possible. Rearranging (6.6) for yields 



= 



V{x-xJ ± \jc1{x-x^)'^ - (V^-c]) (y-y„ 



(6.14) 
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From (6.14) we see that a real solution for t exists only if 
the following inequality holds 



- (V^-Cj) (y-y„)^ a 0. (6.15) 



Hence, a solution exists only within the region defined by 



X a + 



-1 

7 



1/ 2 



(y - y„) for y > y„ , 



(6.16a) 



and 



X 2 x„ + 



~2 



_ - 1 



1/2 



(y« - y) for y < y„ 



(6.16b) 



since T+ and t_ are both positive, we see from (6.7) that t 
must be between and t., ie., (t.stst + ) . Therefore, s must 
be less than or equal to (t-T.). 

To solve Aj(x,y,t) for as t -* », we start with 
(6.3). Integrating (6.3) with respect to X and Y yields 
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A^(x,y, t) = 



2tTC^(V^ C^) -i ( t-T J -s) ( ( t-T.) -s) 



I 



g(s) ds 



(6.17) 



where and t_ are defined in (6.14) . Using (6.1c), we brealc 
(6.17) into the following parts 



AJx,y, t) = , 



(6.18a) 



where 



J ^ 1 r sr(s) ds 

27 TCj(V'-Cj) 1 -s)((t-T.) -s) 



(6 . 18b) 



and 



I. = 






ds 



27 TC^(c^ t-rj -s)(( t-T.) -s) 



(6 . 18c) 



Looking at the behavior of I 3 as t -» 00 , we see from 
(6.18b) that I 3 is 0(l/t). Therefore, we can neglect this 
portion of Aj(x,y,t). The evaluation of I 4 is similar to the 
analysis as developed for Ij: 
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1 



1/2 



= 



27rc,(\^2_^2) 1/2 ^ 



In 1 



( t-s„-a] 



( t-s-a) 



- 1 



-1 


1 1 
In 






' b 


[ J J 



(6.19) 



Here a = (t + + t.)/2 and = (t + -t.)^/4. For long time, I 4 is 
0 (ln t) . 

In the case where , we see that the long time 
solution for Aj(x,y,t) is 0(ln t) . Next, we examine the 
resonant case where = V^. 

B. C,^ = 

Previously, we developed the long time solution for the 
cases when c^^ > and with the interesting result that 
in both cases the solution had growth of 0(ln t) . The reason 
for this growth is the failure of Huygens' principle in two 
dimensions. Though ln(t) tends to infinity for long time, it 
is a function of extremely slow growth, and for the times we 
are interested in, it can be considered as a bounded function. 

In the resonant case, we are concerned with the behavior 
of As(x,y,t) as t ^ 00 for the case when the speed of the 
submarine is equal to the long-wave phase speed (c^ = V^) . 
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To determine if a solution to (6.3) is possible for = 
V^, we must reexamine (6.5). From (6.5) we see that a 
solution exists only if s satisfies the following inequality 



s s t 



(x-xj2 + (y-yj^ 
2V{x-xJ 



for V(x-x^) > 0. 



( 6 . 20 ) 



Integrating (6.3) with respect to X and Y yields 
A^(x,y, t) = 

1 M g{s) ds .21) 

° ^2V(x-Xo) (t-s) - [(x-x„)^ + (y-y„)^] 

Here A = [t - ( (x-x„) (y-y^,) ^) / (2V (x-x„) ) ] . Using (6.1c), we 

brealc (6.21) into the following parts 



A,(x,y, t) = Jj + Ig , 



(6.22a) 



where 
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I. 



1 




g{s) ds 



(6.22b) 



27TC, 



2V'(x-Xo) (t-s) - [(x-x^)^ + (y-y^)^] 



and 



1 




A 



(6.22c) 



2nc. 



Here A = [t - ( (x-x„) (y-y„) ^) / (2V {x-x„) ) ] . 

Looking at the behavior of Ij as t ^ we see from 

(6.22b) that for long time, I 5 is 0(l/t'^^). Therefore we can 
neglect this portion of A 3 (x,y,t) . Integrating directly, we 
obtain 



From (6.23) we see that is 0(t'^^) as t -* 00 . Therefore, we 
see that the solution for As(x,y,t) at resonance grows like 
0 (t‘'^) for long time. 



I. 




(6.23) 



27TC,. 



V(X-X^) 
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To derive the long time solution (t ») for Aj(x,y,t), we 
began by introducing a moving point source for the topography, 
fj(x+Vs,y, s) , with the knowledge that if a solution was 
attainable by using (6.1c) , we could solve (5.23) for any form 
of the topography by the principle of superposition. We first 
looked at the case where c^^ > and determined that the long 
time solution for A,(x,y,t) was 0(ln t) . Next, we looked at 
the case where and saw that a solution existed only 
within a certain region defined by (6.16a,b). The long time 
solution for A 5 (x,y,t) in this case was also 0(ln t) . We 
attribute the 0(ln t) growth to the failure of Huygens' 
principle in two dimensions. For the resonant case where c^^ 
= V^, we found the solution for A 5 (x,y,t) to be 0(t’'^) for long 
time. From this analysis, we see in this problem that if a 
submarine patrols at or near a long -wave phase speed a 
significant resonant could occur for long times. 
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VII. CONCLUSIONS 



This thesis has presented a theoretical study of the flow 
of a stratified fluid over variable topography as a model of 
the resonant interaction of a submarine's wake with a 
stratified fluid. We considered the case in which a submarine 
patrols at a speed near the long -wave phase speed of an 
internal wave. Internal waves are generated by the flow 
around the submarine ' s wake . 

The long time solution obtained in our model for the cases 
when Cj^ > and c^^ < had growth of 0(ln t) . (see (6.13) and 
(6.19)). We attribute the 0(ln t) growth to the failure of 
Huygens' principle in two dimensions. For the resonant case 
where = V^, we found the solution to be 0(t'^^) for long 
time. (see (6.23)). It is important to note that for the 
time scales we are concerned with, both long time solutions 
for the nonresonant cases are bounded and have extremely slow 
growth. Therefore it is the resonant case which is important. 

The growth in time occurring in the resonant case 
indicates that perhaps nonlinear effects, balanced by 
dispersion, could cause significant upstream disturbances. 
From these disturbances, the possibility exists of detecting 
a submarine by the internal waves its wake generates. 
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Further studies in this problem should begin by rescaling 
the asymptotic expansions to balance the leading order 
quadratic nonlinear terms with the topographic forcing. This 
will achieve a balance between nonlinearity and dispersion. 
We expect the solution to this problem to show that 
significant upstream disturbances could occur when the 
submarine patrols at a speed at or near a long-wave phase 
speed. 

If internal waves were generated only by underwater moving 
objects then the problem of detecting the submarine by the 
internal waves its wake generates would be made simpler, but 
this is not the case. Internal waves are a natural occurring 
phenomenon in the oceans. This complicates the problem 
because in order to detect the submarine by the internal waves 
that it generates, the ambient internal wave patterns at the 
location of the submarine must be known. The problem of 
discerning ambient internal waves from internal waves 
generated by the submarine’s wake is similar to the problem of 
detecting a submarine passively (such as with a tail) . The 
sonar operator filters out the ambient noises or background 
noises of the ocean from the incoming signal to detect the 
acoustic signature of the submarine. From further 
understanding of the ocean's naturally occurring internal 
waves, we may one day be able to discern internal waves 
generated by a submarine from those of the ocean. 
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